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Magnetic properties of the spin doughnut M075V20, which is a frustrated magnetic cluster 
with twenty spins, are studied by using the finite-temperature Lanczos method. Our model 
Hamiltonian consists of one ten-spin ring and two five-spin rings with nearest-neighbor anti- 
ferromagnetic couplings J and J". The five-spin rings are coupled to the ten-spin ring by an 
antiferromagnetic exchange interaction J' so as to form isosceles triangles of one J-bond and 
two J'-bonds. Fitting the theoretical magnetic susceptibility with the experimental one mea- 
sured by Miiller et al, the exchange parameters are estimated as J = 388 K, J' = 163 K and 
J" — 81 K. The most interesting feature of this system is that the two five-spin rings behave 
like two free-spins at low temperatures in despite of J' > J". We also calculate temperature 
dependence of the specific heat, and find three peaks at T ~ J, J" and J'/100. 
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A molecular quantum-spin icosidodecahedron was 
found to be realized in a kind of giant molybdenum- 
oxide-based spherical clusters, and has provided a new 
stage of the investigation of frustrated spin systems. 
The term icosidodecahedron means a polygon formed by 
twenty triangles and twelve pentagons. Typical exam- 
ples of molecular quantum-spin icosidodecahedrons are 
Mo72Fe3 U 1 and M072V30, 2 ' 3 which are often called spin 
balls. 

The S = 5/2 spin ball Mo72Fe3o has been studied ex- 
tensively by many authors. 4-9 The main concern of this 
material is to study the crossover between quantum and 
classical regions. In addition, the antiferromagnetic inter- 
action, J ~ 1.6 K, is so small that one can saturate the 
magnetization experimentally by a magnetic field ~ 20 
T. As for the 5 = 1/2 spin ball M072V30, Miiller et al. 
measured the spin susceptibility, which asymptotically 
approaches zero as T — > 0, featuring a maximum at 11.0 
K. 10 They also analyzed this experimental susceptibility 
on the basis of the antiferromagnetic Heisenberg model 
with nearest-neighbor coupling J. 11 By using quantum 
Monte Carlo method to study the high-temperature be- 
havior for T > 120 K, they obtained J = 245K and 
the Lande factor of g = 1.95. However, due to the nega- 
tive sign problem of the method, it is not clear whether 
the lower-temperature susceptibility can be reproduced 
within the nearest-neighbor Heisenberg model or not. 

Replacing ten V 4+ in top and bottom of the spin ball 
M072V30 with neutrality ions, we obtain the S — 1/2 
spin doughnut M075V20, 2 ' 3 which is the main subject 
in this paper. The susceptibility x of the spin doughnut 
was measured by Miiller et al. from T = 300 K to 2 K. A 
quite interesting behavior was found in the experiment: 
x(T) increases as temperature decreases, and shows no 
maximum down to 2 K. Especially, they pointed out that 
the susceptibility follows the Curie law, x — const/T, 



below 20 K. The magnitude of the Curie constant is as 
large as that of two free-spins par spin doughnut. In this 
letter, we assume that the low-temperature Curie law 
is an inherent property of M075 V20, and study how the 
free-spin like behavior emerges. 

Magnetism of M075V20 is expected to be described by 
the following model Hamiltonian: 

L/2 L 

H =J S , 2i-i • S , 2i+i + J' ^2 Si ■ Si+i 

i=l i=l 
L/2 

+J"j2s 2i -S 2i +4 (1) 

with L = 20 (see Fig. 1). As seen in Fig. 1, J- and J"- 
bonds, respectively, construct a ten-spin ring and two 
five-spin rings, and a spin on the five-spin ring is con- 
nected to two spins on the ten-spin ring by J'-bonds. 

This model was proposed for the first time by Miiller 
et el. in their theoretical analysis of the experimental 
susceptibility. 3 Omitting J"-bonds first, they calculated 
temperature dependence of the susceptibility for L = 6- 
16 systems, which are smaller than L = 20 due to the 
limitation of their computer memory, and found that the 
best description of the experimental data was obtained 
for L — 12. It was reported that the low temperature 
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Fig. 1. Schematic representation of the Heisenberg model defined 
in eq. (1) with L = 20 for the spin doughnut Mo75V20- 
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Curie law due to two uncorrelated spins was reproduced 
for a narrow range around J' = 0.55J. They finally in- 
troduced J"-bonds to improve the agreement between 
the theory and experiments, and resulted in J = 288K, 
J' = 0.55 J, J" = 0.20 J and g = 1.95. 

In this letter, we calculate the eigenvalue distribution 
function of eq. (1) with L = 20 to study the suscepti- 
bility and specific heat of the spin doughnut Mo7 5 V 2 o- 
We are going to show that the present model shows the 
low-temperature Curie law for a wide parameter region, 
which is not restricted to the region around J' = 0.55 J. 
The origin of the low-temperature Curie law is that the 
five-spin rings are in the doublet states at temperatures 
below J". As described later, our theoretical analysis 
of the experimental susceptibility of M075V20 leads to 
J' = 0.42J and J" = 0.21J. Note that J' > J" is natu- 
rally expected from possible exchange paths in Mo75"V2o- 
It should be stressed that, in spite of J' > J", there ap- 
pears the doublet behavior of the five-spin rings in ther- 
modynamics. Because of the zigzag form of J'-bonds, 
the five-spin rings are possible to be almost independent 
from the ten-spin ring in MO75V20- 

We here describe our calculation method. The spin 
doughnut has twenty 5=1/2 operators, so it is difficult 
to calculate all eigenvalues by the Householder method. 
An alternating way is to calculate the eigenvalue distri- 
bution function (EvDF) 

P (lo) = Tr 5(u> - H) (2) 

by using the finite-temperature Lanczos method. 12 The 
EvDF can be rewritten as 

M 

^ N M 

= J2J2lm M (n\G(LJ + iO)\n)M, (3) 

11 M n=l 

where {|n)M} denotes the Ising basis with the total S z of 
M, Nm is the dimension of the M-subspace, and Q(z) — 
(z — Ti)^ 1 is Green's function. In the finite-temperature 
Lanczos method, the sum over n is estimated by sam- 
pling random vectors {<()} in the M-subspace, and the 
matrix elements m(4>\G(z)\4>)m are calculated by using 
the continued fraction form with Lanczos coefficients of 
the Hamiltonian matrix. Once the EvDF is obtained, ex- 
pectation values of H l and (S£ ot ) 1 are calculated via 

(n l )=zj_ due-^Jp{w), (4) 
{{SU l ) = - / dwe-^MWM, (5) 

J -°° M 

where Z = dtoe"^ p(u), and thus we obtain the sus- 
ceptibility X = L- 1 [((5 t z ot ) 2 } - (5 t z ot ) 2 ]/T and the specific 
heat C = L-^H 2 ) - {H) 2 ]/T 2 per spin. It has been 
known that this procedure gives fairly good results of 
thermodynamic quantities. 12 ' 13 

Before proceeding to our calculated results, let us 
make a rough estimation of a = J' / J and b = J" /J 
for Mo75"V2o- We expect b < a, because the exchange 



paths of J" are longer than those of J'. So we first con- 
sider the ground state phase diagram of the sawtooth 
lattice Heisenberg model with J and J', omitting J"- 
bonds. We call spins on tips of the sawtooth "tip spins" 
and the others "bottom spins". At a = 1, the ground 
state is the exact dimer ground state, where tip and bot- 
tom spin pairs are in the singlet dimer state. 14 By using 
the level-spcctroscopy method, 15 Blundell and Nuncz- 
Regueiro studied the quantum phase transition point be- 
tween the dimer phase and the spin fluid phase in the re- 
gion of a > 1, and showed that the transition point exists 
at a = 2.052. 16 They also studied the region of a < 1, 
and found that the dimer phase ends at a = 0.653. 16 In 
the ground state for a > 0.653, the dimer or spin-fluid 
state, a tip spin correlates with a bottom spin. Then, 
we expect small b term gives just a perturbation to the 
system, and the low-temperature Curie law is not ob- 
tained. Therefore, in order to get the low-temperature 
Curie law, we concentrate on the region of a <~ 0.6, 
where b terms may play an essential role, not a weak 
perturbation. If we suppose the low-temperature Curie 
law is obtained, the onset temperature should be domi- 
nated by the strength of J"-bonds. As suggested by the 
study of the spin ball, we expect the magnitude of J in 
the spin doughnut is several hundred K. The experiment 
on the spin doughnut has indicated that the onset tem- 
perature of the low-temperature Curie law is about 20 
K, which leads to b ~ 0.1. 

The crucial point of our scenario is whether the ef- 
fect of J"-bonds indeed dominates over that of stronger 
J'-bonds. In order to check this point, we show our calcu- 
lated results of \T per spin for several values of a and b in 
Fig. 2. We find in this figure that the susceptibility shows 
the Curie law of two-uncorrelated spins, %T = 0.025, at 
low temperatures for all parameter sets used in this cal- 
culation. 

We turn to analysis of the experimental susceptibility 
to estimate the exchange parameters for Mo7sV2o- We 
intend to make two types of fittings, which we call "case 
1" and "case 2", respectively. In case 1, we use the ex- 
perimental data of the fresh sample in Fig. 16 of Ref. 2. 
The sample is assumed to contain no magnetic impuri- 
ties, and the Lande factor is determined by the condition 
that the low temperature tail of the experimental data 



_ 1 1 1 1 1 

1 a=0.45 


1 


1 


b=a/2 






■ b=a/3 






' a=0.42 




w - 


_ -- b=al1 




--- b=a/3 






a=0.40 






- - b=a/2 




/ ; 


- — - b=ali 






■ 0=0.35 






- - *=o/2 































Fig. 2. Overall behavior of %T par spin for several values of a = 
J' I J and b = J" /J. A plateau in \T a * low temperatures is seen 
in all data. 
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agrees with the Curie law of two uncorrelated-spins. In 
case 2, we take into account the existence of magnetic 
impurities with 5 = 1/2 and the Lande factor is fixed 
to g = 2. When g = 2 is assumed, the low temperature 
tails of the fresh sample data and the aged sample data 
in Fig. 16 of Ref. 2 are somewhat larger than that of 
the Curie law of two uncorrelated-spins. We use those 
differences to estimate concentrations of the magnetic 
impurities. 

Following the fitting procedure of case 1, we obtain 
g = 2.24 and find the best fit is achieved at J = 575 
K, a = 0.35 (J' = 201 K) and b = 0.21 (J" = 101 K). 
For the spin ball (SB) M072V30, the exchange constant 
Jsb and the Lande factor <?sb have been estimated as 
Jsb = 245 K and gsB = 1-95, respectively. 10 The discrep- 
ancy between those and the present results of J = 575 
K and g = 2.24 for the spin doughnut may be too large. 
In Fig. 3(a), the calculated results of \T and \ are com- 
pared with the experimental result. We find in this figure 
that disagreement is pronounced at intermediate temper- 
atures, 0.04 < T/J < 0.15. 

When we use the procedure of case 2, the impurity 
concentrations of the fresh and aged samples are, respec- 
tively, estimated as 1.3% and 4.0%. We also found that 
both samples give almost the same curve, after subtract- 
ing the impurity contribution from the raw data of the 
susceptibility. Using this susceptibility data, we obtain 
J = 388 K, a = 0.42 (J' = 163 K) and b = 0.21 (J" = 81 
K) via the fitting procedure. The value of J = 388 K is 
closer to Jsb = 245 K than case 1. The theoretical and 
experimental curves are shown in Fig. 3(b). The agree- 
ment is improved compared with case 1. Therefore we 
will use the results of case 2 as our estimation of the 
model parameters for Mo7 5 V2o below. 

We here comment on the magnetic impurities. In the 
experiment of the spin ball M072V30, the fresh sample 

(a) case 1 (b) case 2 
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Fig. 3. Comparison between the calculated and experimental sus- 
ceptibilities per spin. In (a) case 1, we assume there are no mag- 
netic impurities. In (b) case 2, we fix the Lande factor as g = 2.00, 
and use the impurity concentration, ni mp , as a fitting parameter. 
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gives almost the same susceptibility to the aged sam- 
ple and no Curie term due to magnetic impurities is 
observed. 10 On the other hand, for the spin doughnut 
M075V20 the difference of the aged and fresh samples 
is visible in the experiment. 3 As mentioned above, we 
have found that the difference between the two sam- 
ples is resolved by subtracting the impurity contributions 
and better agreement between the theory and experi- 
ment is obtained for the subtracted data than the raw 
data. These facts may suggest that the main difference 
of the aged and fresh samples is in the impurity concen- 
tration, although the present authors have no idea about 
the mechanism of the aging process of Mo75V2o- 

In order to get further insight into the magnetism of 
the spin doughnut, we study the temperature depen- 
dence of specific heat. The calculated results are shown 
in Fig. 4, where we find a three-peak structure. The spe- 
cific heat for T/J > 0.03 is approximately given by the 
sum of those of two five-spin rings and one ten-spin ring, 
and thus the peak temperatures of the highest- and mid- 
temperature peaks scale J and J", respectively. At low 
temperatures well below the mid-temperature peak posi- 
tion, two five-spin rings are expected to be in the doublet 
states, which is the origin of the low-temperature Curie 
law. As seen in Fig. 4, the lowest-temperature peak po- 
sition depends on the value of a, and it exists around 
T = J'/IOO. Below this temperature, the J'-bonds get 
two five-spin rings and one ten-spin ring to be coupled 
each other to form a singlet ground state. 

It is instructive to look into the energy histograms 
Pm{oj) for the total S z of M = and 1. For our model 
with J' = 0.42J and J" = 0.21J, the ground state en- 
ergy is given by E g — — 5.5503J. Low energy parts of 
Pm(^) are shown in Fig. 5(a) for uj ~ E g <^ 2 x J" and 
in Fig. 5(b) for u - E g <~ 2 x J'/IOO, where the energy 
mesh, Auj, of the histogram is chosen as Aui = 10~ 5 J. 
We find in Fig. 5(a) that there is a group of states at 
oj ~ Eg. Only the states in this group contribute to ther- 
modynamics at low temperatures well below J", because 
energies of other excitations are higher than those by 
0.17 J ~ J". As shown in Fig. 5(b), this group consists of 
four singlet and four triplet states. Noting that the five- 
spin ring has two-fold degenerate doublet ground states 
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Fig. 4. Temperature dependence of C(T) of the spin doughnut 
for a = 0.45, 0.42 and 0.45, where b = 0.21 is fixed. We used 
J = 388 K to obtain the temperature in K. 
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Fig. 6. Several low-energy eigenvalues for Stot = and 1, as a 
function of a = J' /J. The inset shows the singlet-triplet gap, 
A s t, versus a 2 . 



Fig. 5. Energy histograms pm(u) with total S z of M = and 
1, where energy regions are chosen as (a) u> <~ E g + 2J" and 



(b) u <r 



2J'/100. In the lower panel, s(Z) (t(Z)) denotes 



singlet (triplet) state with Z-fold degeneracy. 



and the ten-spin ring has a singlet ground state, we inter- 
pret the four singlet and four triplet states as originating 
from direct products of ground states of the five- and 
ten-spin rings. 

We here study the low energy singlet and triplet states 
as a function of a = J' /J. We calculate the lowest two 
eigenvalues in the singlet sector and the lowest three 
eigenvalues in the triplet sector by the standard Lanc- 
zos method, and the results are shown in Fig. 6. As ex- 
pected, we find that the energy eigenvalues reduced to 
the sum of ground state energies of five-spin and ten- 
spin rings as a — ► 0, and the energy splitting does not 
grow so much even for a > b. It is natural to expect that 
this feature is originating from the geometrical structure 
of the J'-bond, by which a spin on five-spin ring is con- 
nected to two spins on ten-spin ring. In order to clarify 
this point, we plot the singlet-triplet gap A st , which is 
defined by energy difference between the lowest eigen- 
values in the singlet and triplet sectors, as a function 
of a 2 in the inset of Fig. 6. This plot gives the asymp- 
totic form of A st /J ~ 0.033a 2 . Thus, the smallness of 
the energy splitting is coming from the fact of the lack 
of the first order term and the smallness of perturba- 
tion coefficients for higher order terms. We have to look 
into perturbation processes to understand the smallness 
of the perturbation coefficients. However, such analysis 
remains as a future problem. 

In summary, we have made a theoretical analysis of 
the susceptibility of the spin doughnut M075V20, where 
the impurity concentration has been estimated as 1.3% 
(4.0%) for the fresh (aged) sample and the exchange pa- 
rameters J = 388 K, J' = 163 K and J" = 81 K have 
been extracted from the experimental data. On the basis 
of various calculations, we have discussed that the spin 
doughnut should be viewed as a weak coupling system of 
two five-spin rings and one ten-spin ring, for which ap- 
pearance of the low-temperature Curie law is quite nat- 
ural because of the doublet states in the five-spin ring. 



We have also predicted that the low-temperature Curie 
law ends around T = 1 K, together with the lowest tem- 
perature peak of the specific heat. It is desired to check 
those predictions experimentally, which judges the va- 
lidity of our theory. It is also interesting to observe the 
low-temperature specific heat under magnetic fields, be- 
cause both the singlet and triplet states contribute the 
lowest temperature peak of the specific heat and the en- 
ergy differences are so small that modest values of the 
magnetic field can get the specific heat to change. 
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